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ABSTRACT 



We present multicolor Hubble Space Telescope (HST) WFPC2 broadband 
observations of the Type Ic SN 19941 obtained ~280 d after maximum light. 
We measure the brightness of the SN and, relying on the detailed spectroscopic 
database of SN 19941, we transform the ground-based photometry obtained at 
early times to the HST photometric system, deriving light curves for the WFPC2 
FA39W, F555W, F675W, and F8UW passbands that extend from 7 d before 
to 280 d after maximum. We use the multicolor photometry to build a quasi- 
bolometric light curve of SN 19941, and compare it with similarly constructed 
light curves of other supernovae. In doing so, we propose and test a scaling in 
energy and time that allows for a more meaningful comparison of the exponen- 
tial tails of different events. Through comparison with models, we find that the 
late-time light curve of SN 19941 is consistent with that of spherically symmetric 
ejecta in homologous expansion, for which the ability to trap the 7-rays pro- 
duced by the radioactive decay of 56 Co diminishes roughly as the inverse of time 
squared. We also find that by the time of the HST photometry, the light curve 
was significantly energized by the annihilation of positrons. 



Subject headings: Supernovae 
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1. Introduction 

Late-time light curves are an important piece of evidence for understanding supernova 
(SN) explosions. Conditioned by the always diminishing optical depth in the ejecta a time 
is always reached when the energy from radioactive decay is radiated on a time scale short 
compared with the dynamical time scale. Radiative equilibrium is established and there is a 
simple connection between energy produced, trapped, and radiated. The simplified physics 
permits a fairly direct connection between the light curve and the dynamical evolution 
of the expanding ejecta. This gives a more robust handle on some of the basic physical 
parameters of the explosion. How fast these late times are reached depend on the mass to 
energy ratio of the ejecta. For supernovae (SNe) without massive H envelopes they start a 
few weeks after maximum light. 

The main physical processes at play are the emission of 7-rays and positrons from 



1969|) 



radioactive decays in the 56 Ni — > 56 Co — > 56 Fe chain (IColgate &: McKed 
their interaction with the ejecta, and the spectrum of the radiation produced by the 
thermalization processes and the radiative transfer in the expanding ejecta. Some knowledge 
of the spectrum is important because ideal bolometric light curves are seldom observed, 
and astrophysicists generally must deal with observations that cover only a fraction of 
the radiated energy. The decay of 56 Ni into 56 Co releases energy in the form of 7-rays, 
with an e-folding time of 8.8 d, while the decay of 56 Co into 56 Fe releases ~ 96% of the 
energy in 7-rays and th e remaining ~ 4% as positrons, with an e-folding time of 111.3 d 



( IDiehl fc Timmes 



1998). Given a density structure and a chemical composition for the 
expanding ejecta, the 7-ray transport problem is determined. Positron interaction with the 
ejecta, on the other hand, strongly d epends on the presence, and geometry, of magnetic 



fields ( jRuiz-Lapuente fc Spruit 



1998|) . These are, and will probably remain, essentially free 



parameters. The escape of positrons from Type la supernovae (SNe la) may be relevant 
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to the diffuse 511 keV positron annihilation radiation from the Galactic bulge region 



( IKnodlseder et al 



20051 ) 



Among SNe that explode thro ugh the gravitati onal collapse of the core, there are two 



distinct scenarios. For SNe II (see 



Filippenko 



19971 . for a review of SN types), associated 



with explosions of massive stars that have retained most of their H /He envelopes, and 
which are endowed with a relatively large mass-to-energy ratio, the ejecta are able first 
to trap the UV and optical-near IR photons for ~100 days, and then to trap the 7-rays 
from radioactive decay s very efficiently for hundreds of days (e.g. 



Elmhamdi et al 



Clocchiatti et al. 



1996a; 



20031 ). Since they trap essentially all of the energy from the radioactive 



decay of 56 Co during the time when both 7-rays and positrons are being produced, their 
light curves do not provide an effective diagnostic for the balance of 7-ray versus positron 
energy deposition. Something similar happens for many (though not necessarily all) Type lb 
and Type lib explosions. 

SNe of Type Ic, on the other hand, associated with the explosion of stars that were 
massive while on the main sequence but that through individual and/or interactive binary 
evolution have lost all of the H and most of, if not all, the He of the outer layers, 
present more favorable conditions to study this problem. Their ejecta have a mass, and/or 
mass-to-energy ratio, small enough to allow a major fraction of the 7-rays from 56 Co to 
escape, potentially providing a cleaner window into the details of positron energy deposition. 

The Type Ic SN 19941 in M51 presented a remarkable opportunity to build a useful 
late-time light curve for ejecta having both a low mass and a low mass-to-energy ratio. It 
exploded in a nearby galaxy and thus reached a bright apparent magnitude at maximum. 
Being in a host galaxy that is routinely observed by amateur astronomers meant, in addition, 
that t he SN was spotted, reported, and spectrqscopically confirmed w ell before maximum 



light (IPuckett et al 



1994 



Schmidt et al. 



1994 : 



Clocchiatti et al. 



19941 ) . Consequently, it 
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received extensive coverage both near maximum light and at later times. 



An optical throu gh near-in f rared (IR) light curve extending ~50 d after maximum 



light was prepared by 



Schmidt 



maximum, was presented by 



(199 41), and a more extensive one, reaching ~120 d after 



Richmond et al. 



(119961 ) . The SN was, in addition, the tar get 



of many other prog rams, including radio and X-ray observations (IRupen et al. 



Immler et al. 



1994a: 



20021 ) . and high- resolution spectroscopy of narrow absorption lines produc ed 



by significant amounts of foreground interstellar matter in M51 (IHo fe Filippe nko 



Ext ensive spectroscop i c data sets were collected and presented by 



and 



Clocchiatti et al. 



(Il996bl ). Finally, the light curve derived by 



'ilippenko et al. 



Schmidt 



1995) 



(1995) 



(jl994J) was 



the subject of detailed theoretical studies to set constraints on the physical parameters of 



the progenitor of SN 19941 and its 



Iwamoto et al 



1994 



Woosley et al. 



evolu tionary path to explosion (jNomoto et al. 



19951 ). 



1994 : 



In this paper, we present Hubble Space Telescope (HST) multi-wavelength broadband 
observations of SN 19941, taken ~280 d after maximum light (~300 days afte r explo sive 



Richmond et al 



(119961 ) . we 



nucleosynthesis). Combining these observations with those of 
build an optical through near-IR light curve that, via detailed modeling, can be used to set 
constraints on the ejected mass and its dynamical and thermal evolution. We introduce the 
new observations and photometry in § [21 In § [3] we describe the problems that need to be 
solved in order to combine the ground-based and space-borne observations, as well as the 
steps we took to construct our optical through near-IR light curve, and in §H]we present the 
results. Our discussion in § [5] gives some physical interpretation based both on simplified 
models of the ejecta and 7-ray radiative transf er, and presents a revision and extension of 



the detailed model C021 (jlwamoto et al. 



19941 ) . We summarize our conclusions in § [61 
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2. Late-Time Observations and Photometry 

The late-time data on SN 19941 were obtained with HST under the Supernova 
INtensive Studies (SINS) General Observer Program. On 1995 Jan. 15 (UT dates are used 
throughout this paper), about 280 d after maximum light, the WFPC2 on board HST 
was used to take five images of the SN in optical and near-IR passbands, totaling 3200 s 
of exposure time. The pointing of the telescope placed the SN on the Planetary Camera 
chip. Table [TJ gives additional details of the observations. The images were reduced using 
the "On-The-Fly Reprocessing" pipeline, provided by the Space Telescope Science Institute 
(STScI). 

Careful inspection of the images revealed the SN at the expected position, amidst the 
complex structure of one of the inner arms of M51. The right ascension and declination 
provided by the world coordinate system of the images was within 1 " of the position fou nd 



1994bJ). 



by high-resolution radio interferometry with the Very Large Array (IRupen et al. 
The identification was confirmed by analysis of two earlier sets of ultraviolet (UV) images 
also obtained with HST. They were obtained with WFPC2, when the SN was close to 
maximum light, using the F336W passband. One image, with an exposure time of 180 s, 
was taken on 1994 April 18, and a set of three images, with a combined exposure time of 
1200 s, was obtained on 1994 May 12. The object identified as SN 19941 appears with very 
good signal in all images, and its flux diminishes by about a magnitude between the two 
epochs. A finder chart showing SN 19941 in the HST UV images is given in Figure [TJ and 
stamps of the region around the SN in the UV, optical, and near-IR images are given in 
Figure El 

The background against which the SN appeared is complex, but the exquisite resolution 
of PCI placed the SN more than 15 pixels away from the closest source. The background 
light in the immediate vicinity of the SN appears to be relatively flat. Point spread function 
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(PSF) fitting photometry of the SN was done us i ng HS TPhot (IDoiphin 



calibrated with the zero points of 



Holtzman et al. 



2000), which is 



(119951 ) by normalizing the PSF magnitude 



to aperture photometry with a 0.5" radius. The adopted version of HSTPhot had zero 
points and charge-transfer efficiency (CTE) correction parameters updated in 2002 October. 
The HSTPhot photometry of SN 19941 is given in Table [H 



3. Merging the Datasets and Estimating Extinction 



If the late-time photometry of S N 19941 presented in th e previous section is combined 



with the ground-based photometry of 



Richmond et al. 



(119961 ). the result will be a multicolor 



light curve extending from ~ 9 d before to ~ 280 d after maximum light. Very few 
SNe Ic have a multicolor light curve extending so long. Combining the datasets requires 
transforming the broadband observations between the different photometric systems used, 
and applying a correction for interstellar extinction. 



The ground based observations of 



Richmond et al. 



(119961 ) were taken with three 



different c ombina t ions o f instrument and telescopes but all of them referred to the standard 



system of iBesselll (1l990l ). We took the data from their Tables 4, 5, and 6, without making 
distinctions between datasets obtained with different telescopes, and selecting the dates for 
which all the B, V, Rc, and Ic bands had been observed. When more than one datapoint 
existed for the same day we combined them by a weighted average. This approach means 
that systematic differences in the calibration of the SN photometry due to the slightly 
different realizations of the standard system are interpreted as an uncertainty in the 
observation. The selected photometry is given in Table El 
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3.1. The Spectral Shape Corrections 

Merging the space-borne and ground-based datasets calls for careful consideration of 
the passbands used. The differences in the SN flux measured by two different realizations of 
the same passban d, evaluated in magnitudes, have been called "spectral shape corre ctions," 



or S-corrections (IStritzinger et al. 



2002 



Krisciunas et al. 



2004 



Pignata et al. 



20041 ). We 



generalize here the use of this name to the difference between the flux measured by a 
ground-based passband and its space-borne counterpart. 

Computing the S-corrections requires precise knowledge of the passbands' throughput 
and, ideally, spectra of SN 19941 at the spe cific dates of the photometry. The spectroscopi c 



evolution of SN 19941 was closely followed (IFilippenko et al 



1995 



Clocchiatti et al. 



1996b) 



although the last spectrum available was taken only ~147 d after maximum light, more than 
130 d before the HST photometry. Regarding the passbands, those of the HST WFPC2 are 
well characterized.], b ut three slightly differe nt realizations of the ground-based BV(RI)c 



system were used by 



Richmond et al 



(119961 ). Hence, as a model for the groun d-based 



system , we computed th e ff-corrections as if th e standard system were that of 



(119901 ) , the one to which iRichmond et al. 



Bcsscll 



(119961 ) attempted to transform. We estimated an 



uncertainty for these S"-corrections by using different realizations of the Bessell passbands 
at assorted optical/near-IR imagers at the Cerro Tololo Interamerican Obs ervatory and the 



Besselll ( I1990l ) 



Kitt Peak National Observatory. Note that we transformed the curves of 
from energy-based to photon-based units, and, for the other ground bassed passbands, we 
factored in a typical CCD quantum efficiency, the spectral response of two aluminum-coated 
surfaces, and the typical response of the sky telluric absorption bands if needed. 

With this database of passbands and spectra, we explored whether it was better 



1 See the URL |ftp: / / ftp.stsci.edu/cdbs| maintained by the STScI. 
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to transform the photometry from the space-borne to the ground-based system, or 
vice-versa. We found that because spectra are available at the times of the ground-based 
observations, it was preferable to convert the ground-based photometry to the space-based 
sy stem. We compu t ed ff- corre ctions for the ground-base d photometry using the spectra 



of 



Filippenko et al. 



119951) and 



Clocchiatti et al. 



(1l996bl ). and passbands as described 



above. We fit smooth polynomials to th eir time variation and u sed them to interpolate the 



S-correction and its uncertainty for the 
are included in Table [2j 



Richmond et al 



(119961 ) observations. The results 



3.2. Interstellar Extinction 



All of the photometric and spectroscopic observations of SN 19941 indicate that it was 



heavily reddened by interstellar matter in M51. Different approaches to estimating /. 



included co mparison with deta i 



and spectra (IBaron et al. 



1996 



ed theoretical models of light curves (jlwamoto et al. 



19991 : 



Millard et al. 



have 



199J) 



19991 1 . and high- resolution spectroscopy 



of the narrow absorptio n lines produced by the gaseous phases of the interstellar matter 



( iHo fc Filippenko 



19951 ). Values of Ay from 0.9 up to 3.1 have been inferr ed, with a 



Richmond et al. 



(119961 ) in adopting 



preferred range roughly between 1 and 2. We follow 
A v = 1.4 ± 0.5 mag. 

This value of Ay must next be translated into extinction for the different broadband 
filters. Typically an interstellar reddening law is adopted and scaled by the value of Ay, 
and the monochromatic fluxes obtained from the photometry are then corrected, applying 
the extinction value that the adopted law prescribes for the effective wavelength of the 
passbands. This approach, however, may lead to systematic errors when applied to objects 
with strong and variable emission lines. The effect is typically small, but it is systematic 
with time, tending to change the overall shape, and/or slope, of a light curve. Since the 
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extinction affecting SN 19941 is large, the issue merits consideration. 



Again, the series of spectra of 



Filippenko et al 



fll995[ ) and 



Clocchiatti et al. 



fll996bh 



allow us to address this issue. Usi ng these spectra, a typical interstellar extinction law 



(ICardelli. Clayton, fc Mathis 



1989), with Ry = 3.1, the transmittance of the passbands 
described earlier, and typical instrumental sensitivities, it is possible to compute the 
extinction affecting the spectra and its variation with time. The estimated extinction values 
and their uncertainties are included in Table [2j 

The lack of a spectrum of SN 19941 some 280 d after maximum affects the estimates 
of interstellar extinction at the time, but the effect is small. Most of the variations that 
strongly affect the S"-corrections at late times (changes of the emission lines in the near-IR) 
have a minor impact on the extinction. 

We computed the extinction for the passbands of the HST WFPC2, fit smooth 
polynomials to the late-time trends, and extrapolated those from the latest observed 
spectrum out to the latest photometric point. To account for the uncertainty of the 
extrapolation, we increased the typical uncertainty of the fits by a factor of three for the 
last point. The extinction computed for the HS T observation is included in Table [TJ 

The uncertainties due to the unknown details of the ground-based standard system 
realization and the lack of a very late-time spectrum are much smaller than the uncer tainty 
of the overall extinction adopted (AAy ps 0.5 mag according to lRichmond et al.lll996l ). The 



changes in the light curves associated with the variation of the broadband extinction with 
time fall well within AAy as well. However, computing the extinction from the evolving 
spectra and the expected passbands is the correct procedure for nonstellar spectra, and 
provides light curves that, for a given velue of Ay, are free from the systematic effects 
related to the time variation of the extinction in each passband. 
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Two additional sources of uncertainty should me mentioned. First, the interstellar 
extinction law may not be the standard one. What can be said in this case is 
that, if the interstellar extinction towards SN 19941 follows the functional form of 



Cardelli. Clayton, fc Mathid (119891 ) with a different value of Ry the difference will be 
smaller than the uncertainty in Ay, even for the more extreme values of Ry that they 
consider. Second, it is possible that at least a fraction of the extinction is produced in the 
local environment of the supernova. Since this local environment will be eventually affected 
by the evolution of the ejecta, a change of the foreground extinction with time would be 
conceivable. With the data available, it is not possible to set quantitative limits to this 
possibility. 



4. The Light Curves 

The resulting light curves of SN 19941 in the space-based system are shown in Figure [31 
while the evolution of the photometric colors is given in Figure HI The late-time point 
is consistent with the decrease in brightness and a flattening of the light curve expected 
from models that fit the post maximum decline (see below), and as suggested by the 
ground-based data earlier than ~80 d after maximum. The (RI)c magnitudes of the last 
two ground-based observations (at ~115 and ~140 d after maximum), and the B and V 
magnitude of the last ground-based observation, are not consistent with this trend. They 
appear to be slightly overluminous. The observed colors indicate little or no change of the 
overall spectral shape since the latest ground-based photometry. 

With all the photometry transformed into the same system, it is possible to convert 
the observed single-bandpass light curves into monochromatic fluxes at the effective 
wavelengths, integrate them over wavelength, and obtain a new version of the lL BV(RI)c" 
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light curve of S IN 



Holtzman et al. 



1994o We used the zero points for the WFPC2 passbands given by 



(119951 ). and integrated the monochromatic fluxes using a simple trapezoid 
rule. We did not make any assumptions regarding the behavior of the flux beyond the 
limiting wavelengths of our coverage; flux beyond the limits of the B and Ic passbands 



was ignored. The distance of M51 was assumed to 



function of planetary nebulae (IFeldmeier et al. 



3e 8.39 Mpc, based on the luminosity 



1993). 



The BV(RI)c light curve of SN 19941 is shown in Figure El toget 



constructed 



Patat et al- 



ight curves of the Type Ic SNe 1990B (IClocchiatti et al 



2001 



Sollerman et al 



2002 



Clocchiatti et al 



rer w ith similarly 



20011), 1998bw_ 



20071 ). and 2002ap flTomita et al. 



2006 ). The latter two have been modele d with ejecta of low mass-to-energy ratio 



( jlwamoto et al. 



1998 



Mazzali et al. 



20021 ). but fairly large ejected mass, and have also 
sometimes been called "hypernovae." As a benchmark of absolute brightness and evolution 
timescales, the BV(RI) C light curve of the normal Type la SN 1992A (N. Suntzeff, 2006, 



private communication) is also given. 



5. Discussion 

5.1. Simple Comparisons 

The light curves in Fig. [5] illustrate the wide range that SNe Ic span in absolute 
brightness, and show that the diversity cuts across observational borders proposed to 
classify them, like the supernova or hypernova categories. Fig. also shows that SN 19941 
was fairly luminous reaching more than 50% of the brightness of normal Type la SN. A 



2 Even though it was calibrated using the space-borne WFPC2 passbands, we prefer to 
name the light curve using the more compact designation of their equivalent ground-based 
passbands. 
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foreground extinction a few tenths of a magnitude larger, which is consistent with model 
light curves, spectra, or high-resolution spectra of the narrow Na D absorption lines, would 
make it as bright as SN 1992A. The BV(RI)c light curves also confirm that, similar to the 
behavior close to maximum light, the time scale for brightness variation of SN 19941 at 
late times was fast. The early exponential tail of SN 19941 is faster than that of the other 
Type Ic SNe, and comparable to that of the Type la SN. 

Qualitative comparison with the late-time light curves of SNe 1990B, 1998bw, and 
2002ap indicates a mismatch between the trend suggested by the last ground-based 
photometric points and the HST observation. SN 19941 shows a sharp downward change 
of slope in the exponential tail after ~140 days, not seen in the other SNe Ic. The rate 
of decay measured from the data between ~80 and ~ 140 d after explosion is 43% slower 
than that measured using the points between ~40 and ~80 d, and the change happens very 
fast, between 80 and 115 d after explosion. By contrast, the decline slope of SN 1990B is 
essentially constant at these times, that of SN 1998bw changes by less than ~28% between 
earlier than 80 d after explosion and later than 250 d after explosion, and that of SN 2002ap 
changes by ~7% in this same time period. 

If the trend suggested by the last two ground-based photometric points were correct, 
then something fairly dramatic happens to the BV(RI)c light curve between ~140 and 
~300 d after maximum light. One possibility is the so-called "inf rared catastrophe" tha t 



shifts energy from the optical/ near-IR passbands into the far-IR (IFransson fe Chevalier 



1989 



Kozma fc Fransson 



19981 ). though it is not expected at such relatively early epochs in 
SNe Ic. Another possibility is that the brightness at 300 days after maximum is affected 
by a different amount of foreground local extinction, as it would happen if the ejecta of 
SN 19941 formed dust. This is, however, not consistent with the small, or nill, change in 
broadband colors. 



-16- 



Another possibility is that the latest ground ba sed photometric points are affected 



Richmond et al. 



(|l996l ) did image 



by contamination from background light from M51. 
matching and subtraction before performing photometry, but the background of SN 19941 
was particularly bright and had a steep gradient over the spatial scales sampled by the 
PSF of the ground-based images. Thus, residual contamination in the ground-based images 
could be present. 



5.2. Comparison with a Simple Model 



When 7-ray radiative tr ansfer can be approximated by a pure absorption process 



({Sutherland fc Wheeler 



19841 ) . the radioactive nickel assumed to be located at the center 
of the ejecta, accounting for a small fraction of the total mass (i.e. a point source model 
for the 7-ray emission), and the SN ejecta modeled by a spherically symmetric distribution 
of matter in free expansion, the optical depth that extinguish es the 7-rays decreases with 



the inverse of time squared (e.g. 



Clocchiatti fc Wheeler 



19971 ). In this case, the energy 



deposited by 7-rays can be approximated by 

£ 7 oc exp (-t/r Ni ) {l - exp [-(T /i) 2 ]} , (1) 

where tn; is the e-folding time for the decay of 56 Co , and To is the characteristic time for 
a given ejecta structure, which depends on its mass M, energy E, density profile, and the 
absorption cross section per unit mass («) that matter presents to 7-rays. For a given 
density profile, To oc (kM 2 / E) 1 / 2 . Soon after maximum light, when radiative equilibrium is 
established, the above assumptions are expected to hold. Then, in addition, between 100% 
and ~96% of the radioactive energy is released in 7-rays, and most of them are trapped 
in the ejecta. In these conditions, E^ gives a good representation of the energy available 
to heat the ejecta and equation [1] can be fitted to bolometric light curves to estimate the 
parameter T . At later times, as the fraction of 7-rays that escape from the ejecta increases, 
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the ~4% of energy released in the form of positrons needs to be considered. In what follows, 
we will consider the two extreme possibilities for thermalization of the positron annihilation 
energy: either they are completely trapped in the ejecta, or they escape. 

Assuming that the BV(RI)c light curve of SN 19941, between ~40 and ~80 d after 
explosion, is a good representation of the bolometric light curve, we fit the light curve to 
determine To ~ 65 d, if the energy from the 56 Co positrons is completely deposited, or 
T ~ 67 d, if the positrons escape from the ejecta without interacting. The BV(RI)c 
luminosity with these fitted simple models is shown in Figure El 

The simple models are very good at fitting the data between ~40 and ~80 d after 
explosion, suggesting some connection with reality, and they do a reasonably good job 
at ~300 d. Since the (bolometric) model implies ~50% more energy than the BV(RI)c 
observation, the data taken at face value would be consistent either with full trapping of 
positrons from 56 Co and a small bolometric correction at ~300 d, or with a smaller (even 
null) bolometric correction and a more complex process of positron transfer that traps 
~50% of their total energy. It should be noted, though, that the simple model without 
positron trapping falls below the observed light curve. The simplest interpretation of this 
observation is that SN 19941 did trap a significant fraction of the positrons. 

The last two points from ground-based photometry (between ~100 and ~120 days 
after explosion) are missed by both fits. We note, in addition, that it is not possible to 
obtain a good fit for the simple models of eq. Q] when all the ground-based observations 
of the exponential tail (between ~40 and ~120 d after explosion) are included, even if 
the HST point at 300 d is excluded. The problem is the fast change in the slope of the 
BV(RT)c light curve between ~80 and ~115 d after explosion. Fitting this break requires 
including an additional free parameter, like two different values of To, which could represent 
a model of two different mass components each of them with different expansion properties. 
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As a comparison, it is possible to successfully fit eq. [T] to the exponential decay of the light 
curves of all the other SNe Ic shown in Fig. 



5.3. Simple Comparison with Hindsight From the Simple Model 

In §5.11 we compared the derived light curves of several SNe. Although it is illustrative 
to see how different SN light curves behave in a plot like that of Fig. [H it is not completely 
correct to compare the evolution of different ejecta on a calendar-time basis. What we 
would like to do is to plot the light curves using units that are more physically meaningful. 

From a few tens to a few hundreds of days after explosive nucleosynthesis, the evolution 
of the bolometric light curves is governed by two main time scales: rc , the e- folding 
time for the radioactive decay of 56 Co into 56 Fe , and To, the characteristic time for the 
decrease of the 7-ray optical depth (eq. [[]). They are, respectively, the time scale for energy 
production and the time scale for energy deposition. This suggests that the natural energy 
scale to compare late-time light curves of different SNe is L Bo i(t n )/ exp (— t n /r Co ), where t n 
is measured since nucleosynthesis, and that the natural scale to compare the time evolution 
of the ejecta is t e /To, where t e should be measured from the start of free expansion. Since, 
at the times of interest, the difference between t n and t e is negligible, either of them may 
be used. We fitted the models of eq. [I] to the light curves of all SNe Ic in Fig. [51 and then 
normalized the energy and time axes according to the previous prescription. The result is 
shown in Figure [71 

It is reassuring to confirm that such a simple scaling brings some "order" to Fig. [5j 
When the energy is normalized by the input 56 Co source, and time by the natural scale for 



3 We note also that eq. [H fails to provide a good fit to the late-time light curve of the 
Type la SN 1992A. 
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7-ray deposition (responsible for most of the energy available to the ejecta from soon after 
maximum until the supply of 56 Co is exhausted), all of the late-time light curves follow 
approximately the same track. It is possible to bring them into close agreement just by 
normalizing the total amount of 56 Ni ejected by the explosion (i.e. the additive constant in 
the logarithm of the energy used in Fig. [5]). 

This comparison shows that, in units of the natural time scale of evolution of the 
7-ray optical depth, the late-time point of SN 19941 presented in this paper (~ 280 d 
after maximum) is even more extreme than those observed for SN 1998bw and SN 2002ap, 
which had light curves extending beyond ~500 d after maximum. Also, it shows that, in 
comparison with Eq. [TJ the three very late light curves are consistent with a flattening after 
~ 3To, a time at which the 7-ray optical depth becomes negligible, and most of the 7-rays 
escape into space. For SN 19941, the exact time at which the flattening occurs is unknown 
due to the lack of observations. However, both at 3T ~ 200 d, and at 300 d corresponding 
to the last observation, a modest fraction of 56 Co remains in the ejecta, and the ~4% of the 
radiated energy that its decay releases in the positron channel becomes a significant source. 

Regarding the late-time ground-based photometry, the scaling of Fig. [7] confirms the 
arguments given before. The last two points depart from the common locus, and probably 
overestimate the luminosity of SN 19941 due to background contamination. 



5.4. Insights From a Detailed Model 



Iwamoto et al. 



(119941 ) presented a model of SN 19941 which, without fine tuning of the 
structure, provides a good fit to the early part of the UBV(RI)c light curve, from ~ 6 d 
before up to ~ 40 d after maximum light. Their preferred configuration is a C+O core 



that would result from the evolution of a massive star in an interacting binary system that 
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underwent two mass-transfer episodes. Depending on the initial parameters of the binary, 



at the moment of explosion the pair consists o 
dwarf, or low-mass main-sequence companion 



the presupernova a nd a neutron star, white 



Nomoto et al. 



199J) 



Iwamoto et al. 



Schmidt 



(1994) 



(119941 ) fitted the light curves of SN 19941 provided by 
with the explosion of a star that entered the main sequence with ~ 15M , became a 
~ 4.OM He star after the first mass transfer, and a 2.1 M C+O core after the second 
(model C021). A simultaneous fit of the monochromatic and UBV(RI)c curves required an 
interstellar reddening of E(B — V) = 0.45 mag, and a distance of 6.92 Mpc (m — M = 29.2 
mag), in good agreement with the values used to determine the BV(RI)c light curve 
presented here. 

Would model C021 fit the very late-time HST photometry? The question is relevant 
because previous attempts at fitting SN light curves over long time spans have uncovered a 
problem. When the models are good at fitting the maximum light they tend to be too dim 
on the tail; conversely, if they are good at fitting the late-time decay, they tend to give a 



maximum broader and/or brighter than observed. Solving this prob 



introducing free paramet ers in the density structure of the model (jEnsman fc Woosley 



1988 



Tomita et al. 



200J). 



em has typically me ant 



Based on the explosion model C021 (jlwamoto et al. 



19941 ). and using the same radiative 



transfer code, we calculated the monochromatic light curves, but here we extend the 
computations to 300 days after explo sion. The assumptio ns an d approximations used in the 



code are described in more detail in (jlwamoto et al. 



19941 ) and (IHoflich. Muller fe Khokhlov 



19931 ). A hypothesis relevant to our discussion below, is that all positrons are thermalized 
locally and cannot escape. Because of the limitations of the code in the treatment of 
forbidden lines and their time dependence, we assumed the energy distribution to be 
unchanged after day 120. This assumption is consistent with the observations because the 
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color indices scarcely change after about 100 d (Tableland Fig. H]). The assumption 
prevents us from testing overall migration of flux toward the IR, as may result from the 
formation of dust. The evolution of the colors, however, suggest that neither dust formation, 
nor other change in the foreground extinction by matter in the circumstellar region, do 
occur in this SN. 

From the monochromatic fluxes, we calculated the magnitudes in each relevant HST 
pass band. We check e d tha t the synthetic F439W and F555W values agree with the B and 



V of llwamoto et al 



()1994j ) within a few hundredths of a magnitude. The small differences 
can be attributed to adjustments in the frequency grid and differences in the HST filter 
functions and the Bessell system, since S'-corrections were not applied in this test. 

From the monochromatic fluxes at the effective wavelength of the passbands, we built 
a BV(RI)c light curve using the same procedure as in the observations (see §1]). We 
compare the theoretical light curves of model C021 with the observations in Fig. [6j Both 
the bolometric and BV(RI)c luminosities are plotted. They are very similar, indicating 
that the model concentrates most of the energy within the passbands observed. Overall, 
the model luminosities agree well with the observations, especially in the post-maximum 
decay, although some details merit further comment. 

At maximum, the model results dimmer than the observations, while the same model 



fits better the UBV(RI)c in llwamoto et al 



(119941 ) . One reason for that is the distance. 



The o 



jserved BV(Rf 



Mpc ( Feldmeier et al. 



r, lig ht curve presented here is based on an observed distance of 8.39 



19971 ). while, for the same foregroun d extinction, C021 is consistent 



with a distance 1.47 Mpc shorter in 



Iwamoto et al. 



(11994I ). The shorter distance will help 



our model BV(RI)c light curve to fit the maximum better, but will give a brighter than 
required post maximum decay. This highlights the importance of using extensive datasets 
when comparing models to observations. 
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The BV(RI)c light curve of C021 shown in Fig. [6] does fit most of the observed 
BV(RI)c light curve of SN 19941. First, the distance independent width of the observed 
maximum, and then the distance dependent brightness of the post-maximum decline 
and exponential decay. Fitting the peak-to-tail contrast will require some additional fine 
tunning. 

A possible hint is that the model UBV(RI)c light curve does a better job at fitting 
the peak-to-tail contrast of the BV(RI)c observations, indicating that the model spectrum 
overestimates the fraction of total flux emerging in the ultraviolet at early times. In any 
case, the problem of C021 is opposite to the traditional peak-to-tail contrast problem: if 
the maximum of the BV(RI)c model is arbitrarily shifted to match the observed one, 
then the model light curve results brighter than the observed, not dimmer. Studying the 
peak-to-tail contrast problem in SN light curves is beyond the scope of this paper, so we 
will let the issue rest there. 

The late-time tail of the detailed model, which includes a realistic approximation to 
the 7-ray transfer and deposition, follows closely the simple deposition function given by 
eq. [I| with full trapping of the positrons from the decay of 56 Co . Since C021 also assumes 
full trapping of positrons, this is an indication that the 7-ray transfer is well represented by 
the simple approximations that went into eq. [U Hence, the conclusion stated above that 
about half of the positrons are able to deposit their energy within the ejecta is supported 
by the more elaborate model as well. 

6. Conclusions 

Through the HST GO SINS program, we observed SN 1994 with WFPC2 when the 
SN was about 280 d past maximum light. We did PSF photometry on the resulting images 
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using HSTPhot flDolphinll2000h . 



We collected published ground-based photometry of SN 19941, applied a reddening 
correction, and transformed the results to the space-based HST WFPC2 photometric 
system. Both the reddening correction and the photometric tr ansformation were bas ed 



on the extensive spectr oscopic database available for the SN ( jFilippenko et al. 



1995 



Clocchiatti et al. 



1996H 1. The resulting light curves in the F439W, F555W, F675W, and 



F81AW passbands span from ~7 d before to ~280 d after maximum light. 



Using the zero points and effective wavelengths of iHoltzman et al.l (Il995l ). we 
transformed the multicolor photometry to monochromatic flux, and then integrated the flux 
over the available wavelength range to construct a u BV(RI)c" light curve. We compared 
the BV(RI)c light curve of SN 19941 with those of other SNe Ic, and, based on a simplified 
model of the ejecta, proposed a scaling in energy and time that provides a more insightful 
comparison of the exponential tails of different SNe. The light curves that reach times larger 
than ~ 3To, where To is the characteristic time for the decay of the 7-ray optical depth, 
show a clear flattening. For rapidly evolving SNe, like SN 19941, a simple interpretation of 
this flattening is that the small percentage of energy that the decay of 56 Co releases in the 
form of positrons becomes a relevant thermal energy source for the ejecta. 



Model C021, originally presented by llwamoto et al.l (119941 ). does a reasonable job at 
fitting both the BV(RI)c and monochromatic light curves at early times. We extended 
the original light-curve calculations up to day 300. We note that the absolute BV(RI)c 
luminosity agrees well with the observations of SN 19941 up to about 80 d after maximum 
light, but overestimates the flux by ~ 50% at day 300. 

Models with full positron trapping indicate that some energy was missed by the 
observations, but models without positron trapping indicate an observed excess of energy. 
The missing (or excess) energy amounts to ~50% of that produced by the annihilation 
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of 56 Co positrons. Although models with no positron trapping and a migration of 
energy toward the mid-to-far IR (e.g., by dust formation) cannot be ruled out, a simpler 
intepretation of these late-time observations is that the ejecta of SN 19941 trapped ~50% 
of the positrons roughly 300 d after the explosion. 

Our interpretation of the BV(RI)c light curve is consistent with the SN evolving as 
a freely expanding envelope whose optical depth to 7-rays decreases, very approximately, 
with the inverse of time squared. Simple models that treat the deposition of 7-rays as a 
pure absorption process, and the more detailed model C021, agree on this. S N 19941 does 
not ne ed a model with two dynamically distinct components, like the one that 



Tomita et al. 



(12006h propose for SN 2002ap. 



We conclude that our very late-time BV(RI)c photometry represents most of the 
energy that the SN produced at this time, implying that the late-time bolometric correction 
of this rapidly evolving SN Ic was not large. This, in turn, would mean that the "infrared 
catastrophe" in this SN does not happen until phases later than ~300 d. 

Although model C021 was not finely tuned to closely match the observed light curves 
of SN 19941, it helps to show that the late-time photometry presented here is extreme in 
terms of the 7-ray optical depth, and agrees with the simple models of eq. [T] in that the 
dominant source of late-time energy originates in positrons. A model built specifically to 
represent this SN could be used to establish a meaningful constraint on the fraction of 
positrons actually trapped in the ejecta, and whether (and when) there is a migration of 
photons from the optical and near-IR to the IR region of the spectrum. 
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Fig. 1. — HST image of SN 19941 and its surroundings. The SN is at the center of the frame, 
inside the circle. North is at the top and East is to the left. The width of the field displayed 
is about 28". The image shown is the combination of the three images taken on 1994 May 
12, with the F336W passband. 
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Fig. 2. — Stamps of the region around SN 19941 in the HST images. North is at the top, 
and East is to the left. The width of each stamp is 4.6". The passbands used are indicated 
by labels. SN 19941 is at the center of each stamp. The top-left stamp corresponds to the 
F336W exposure taken on 1994 April 18, while the top-center stamp is an expanded view 
of the 1994 May 12 image shown in Figure [TJ The dimming of SN 19941 in the 24 elapsed 
days is clear. The other images were taken on 1995 Jan. 15 (see Tabled]). 



31 



i — i — i — I — i — i — i — | — r 

LD O ID 

^ O 



o 

* M - O 

CO 





LD 


LD 


CO 






LD 


CO 


00 




LD 














■ 


• 





o 
o 
cv 



CO 
CO 



H- M » h«H 
H- N WW 



# H * ^ 

/ / / ^ 
^V 1 ••• h^ 1 



H 



O 



o 
Q 

o 
o cd 



o 



o 

(02 



9pr4iu§i3^ SOMA 



Fig. 3.— iJST WFPC2 broadband unreddened photometry of SN 19941 (t 



re last point, at 



JD — 2,449,433 ~ 300), and ground-based B, V, Rc, and Ic photometry ([Richmond et al. 



19961 ) unreddened and translated to the WFPC2 passbands. The vertical marks aboy e the 



lower horizontal axis indi cate the times at which the spectra of 



Clocchiatti et al. 



Filippenko et al. 



(I1995h and 



()1996bl ) were t aken. The shi 



(JD = 2,449,433) estimated by llwamoto et al. 



t in ti me corresponds to the explosion date 



dl99J). 
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Fig. 4. — HST WFPC2 broadband unreddened colors of SN 19941 in the HST passbands 
computed from the photometry of Fi gured The shift in ti me corresponds to the explosion 
date (JD = 2,449,433) estimated bv llwamoto et all |l99J). 
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Fig. 5.— BV{RI) C light curves of the Type Ic SNe 19941, 1990B, 1998bw, and 2002ap, as 
well as of the Type la SN 1992A. The light curve of SN 1990B has been displaced downward 



by 0.5 vertical u nits for clarity- 



fit for SN 1994 1 ( 



( jMazzali et al. 



Iwamoto et al. 



Exp 



osion-date estimates are from a theoretical light-curve 



1994 ). theoretical spectra and light-curve fits for SN 2002ap 



20021 ). and the date of GRB 98 0425 



explosion date was taken from 



Clocchiatti et al 



or SN 1998bw. For SN 1990B the 



(120011 ). who estimated it from comparison 



with SN 1993J, and for SN 1992A the explosion date was assumed to be ~20 d before 
maximum light in B. The slope for full trapping of the 56 Co decay products is also shown. 
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Fig. 6. — BV(RI)c light curve of SN 19941, the simple models given by eq. [fl and model 
C021. 
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Fig. 7.— BV(RI) C late-time light curves of SNe Ic 19941, 1998bw, 2002ap, and 1990B, 
normalized according to the prescriptions of § 15.31 Points mark the observations of SN 19941. 
Solid symbols indicate the region used to fit the time-scale parameter T . A constant was 
added to the light curves of SNe 19941, 1998bw, 2002ap, and 1990B, to make them match 
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Table 1. Late-Time HST Photometry of SN 19941 



HST Dataset 


JD a 


Exp. b 


Passband 


HST Mag. 


Extinction 


U2LQ0201T/2T 


733.21 


1400 


F439W 


21.63 ± 0.02 


1.89 ± 0.033 


U2LQ0203T 


733.22 


600 


F555W 


21.03 ± 0.01 


1.45 ± 0.030 


U2LQ0204T 


733.25 


600 


FQ75W 


20.21 ± 0.01 


1.12 ± 0.010 


U2LQ0205T 


733.26 


600 


F8UW 


19.69 ± 0.01 


0.87 ± 0.010 



a Mean Julian Date of total exposure minus 2,449,000. For phase, subtract 
an additional 433.00. 

b Total exposure time in the passband (in seconds). Two F439W exposures 
of 700 s each were combined in one image. 

c Interstellar extinction in the host galaxy, in the HST passband. 



Table 2. Ground-Based Photometry of SN 19941 



JD a B b V h R c h I c h A B C A v c A Rc c A Ic c SC B d SC v d SC Rc d SC Ic d 



445 


14.41 


± 


0.03 


14.04 


± 


0.02 


13.87 


± 


0.02 


13.77 


± 


0.01 


1 


.83 


1 


.40 


1. 


17 


0.86 


-0.03 


-0.04 


0.05 


-0.01 


447 


13.97 


± 


0.08 


13.39 


± 


0.02 


13.21 


± 


0.02 


13.09 


± 


0.02 


1 


.83 


1 


.40 


1. 


17 


0.86 


-0.04 


-0.04 


0.06 


-0.01 


448 


13.95 


± 


0.02 


13.20 


± 


0.01 


13.02 


± 


0.01 


12.88 


± 


0.01 


1 


.83 


1 


.40 


1. 


17 


0.86 


-0.05 


-0.04 


0.06 


-0.01 


450 


13.68 


± 


0.21 


12.99 


± 


0.02 


12.72 


± 


0.02 


12.60 


± 


0.02 


1 


.83 


1 


.40 


1. 


17 


0.86 


-0.06 


-0.04 


0.06 


0.00 


453 


13.81 


± 


0.07 


12.91 


± 


0.02 


12.62 


± 


0.02 


12.49 


± 


0.02 


1 


.83 


1 


.40 


1. 


17 


0.86 


-0.07 


-0.04 


0.06 


0.00 


454 


14.05 


± 


0.02 


13.02 


± 


0.01 


12.72 


± 


0.01 


12.52 


± 


0.01 


1 


.83 


1 


.40 


1. 


17 


0.86 


-0.07 


-0.04 


0.06 


0.0^ 


455 


14.15 


± 


0.02 


13.11 


± 


0.01 


12.76 


± 


0.01 


12.53 


± 


0.01 


1 


.83 


1 


.40 


1. 


17 


0.86 


-0.08 


-0.04 


0.06 


I 

0.00 


456 


14.44 


± 


0.01 


13.23 


± 


0.01 


12.86 


± 


0.01 


12.61 


± 


0.01 


1 


.83 


1 


.40 


1. 


17 


0.86 


-0.08 


-0.04 


0.06 


0.00 


457 


14.56 


± 


0.02 


13.34 


± 


0.01 


12.94 


± 


0.02 


12.64 


± 


0.01 


1 


.83 


1 


.40 


1. 


17 


0.86 


-0.08 


-0.04 


0.06 


0.00 


458 


14.75 


± 


0.02 


13.46 


± 


0.02 


13.04 


± 


0.02 


12.72 


± 


0.02 


1 


.83 


1 


.41 


1. 


17 


0.86 


-0.08 


-0.04 


0.06 


0.00 


461 


15.28 


± 


0.03 


13.93 


± 


0.02 


13.41 


± 


0.02 


13.00 


± 


0.02 


1 


.83 


1 


.41 


1. 


17 


0.86 


-0.09 


-0.04 


0.06 


0.00 


462 


15.40 


± 


0.04 


14.08 


± 


0.02 


13.53 


± 


0.03 


13.10 


± 


0.02 


1 


.83 


1 


.41 


1. 


17 


0.85 


-0.09 


-0.03 


0.06 


0.00 


470 


16.00 


± 


0.32 


14.92 


± 


0.02 


14.37 


± 


0.03 


13.79 


± 


0.01 


1 


.84 


1 


.41 


1. 


17 


0.85 


-0.09 


-0.03 


0.07 


0.01 


471 


15.97 


± 


0.30 


15.11 


± 


0.02 


14.50 


± 


0.02 


13.96 


± 


0.02 


1 


.84 


1 


.41 


1. 


17 


0.85 


-0.09 


-0.03 


0.07 


0.01 


476 


16.28 


± 


0.10 


15.41 


± 


0.05 


14.80 


± 


0.07 


14.12 


± 


0.01 


1 


.84 


1 


.41 


1. 


17 


0.85 


-0.08 


-0.02 


0.07 


0.01 



Table 2 — Continued 



JD a B h V h R c h Ic h A B C A v c A Rc c A Ic c SC B d SC v d SC Rc d SC Ic d 



482 


16.61 


± 


0.08 


15.56 


± 


0.02 


15.03 


± 


0.03 


14.37 


± 


0.02 


1 


.84 


1 


.41 


1. 


17 


0.85 


-0.06 


-0.02 


0.07 


0.01 


486 


16.56 


± 


0.11 


15.71 


± 


0.03 


15.20 


± 


0.04 


14.53 


± 


0.01 


1 


.84 


1 


.41 


1 


17 


0.85 


-0.05 


-0.01 


0.08 


0.01 


487 


16.78 


± 


0.04 


15.76 


± 


0.02 


15.24 


± 


0.03 


14.54 


± 


0.02 


1 


.84 


1 


.41 


1. 


17 


0.85 


-0.05 


-0.01 


0.08 


0.01 


491 


16.58 


± 


0.16 


15.90 


± 


0.03 


15.41 


± 


0.07 


14.62 


± 


0.02 


1 


.84 


1 


.41 


1. 


17 


0.85 


-0.04 


-0.01 


0.08 


0.01 


494 


16.88 


± 


0.07 


15.94 


± 


0.02 


15.45 


± 


0.03 


14.75 


± 


0.02 


1 


.84 


1 


.41 


1. 


17 


0.85 


-0.03 


-0.01 


0.08 


0.02 


495 


16.86 


± 


0.09 


15.89 


± 


0.02 


15.49 


± 


0.02 


14.80 


± 


0.03 


1 


.84 


1 


.41 


1. 


17 


0.85 


-0.03 


-0.01 


0.08 


0.0^ 


502 


17.21 


± 


0.14 


16.18 


± 


0.07 


15.73 


± 


0.09 


15.00 


± 


0.03 


1 


.84 


1 


.41 


1. 


17 


0.85 


-0.02 


0.00 


0.08 


0.02 


503 


17.06 


± 


0.06 


16.24 


± 


0.02 


15.75 


± 


0.04 


15.04 


± 


0.02 


1 


.84 


1 


.41 


1. 


17 


0.85 


-0.02 


0.00 


0.08 


0.02 


505 


16.93 


± 


0.17 


16.21 


± 


0.04 


15.74 


± 


0.06 


15.08 


± 


0.02 


1 


.84 


1 


.41 


1. 


17 


0.85 


-0.01 


0.00 


0.09 


0.02 


512 


17.21 


± 


0.18 


16.48 


± 


0.04 


15.99 


± 


0.06 


15.30 


± 


0.02 


1 


.84 


1 


.42 


1. 


17 


0.85 


-0.01 


0.01 


0.09 


0.02 


515 


17.48 


± 


0.08 


16.58 


± 


0.04 


15.97 


± 


0.03 


15.37 


± 


0.05 


1 


.84 


1 


.42 


1, 


17 


0.85 


0.00 


0.01 


0.09 


0.02 


548 


18.05 


± 


0.13 


17.23 


± 


0.08 


16.46 


± 


0.08 


15.95 


± 


0.05 


1 


.85 


1 


.42 


1, 


17 


0.86 


0.00 


0.03 


0.11 


0.03 


572 


18.10 


± 


0.18 


17.60 


± 


0.06 


16.93 


± 


0.09 


16.34 


± 


0.05 


1 


.85 


1 


.42 


1. 


17 


0.86 


0.00 


0.04 


0.12 


0.04 



1 Mean Julian Date minus 2,449,000. For phase, subtract an additional 433.00. 



b Photometry from Tables 4, 5, or 6 of 



Richmond et al. 



(I1996h . 



c Interstellar extinction in the host galaxy, in the given ground-based passband. Mean uncertainties (mag) are 5A B = 0.060, 
5A V = 0.010, 5A Rc = 0.003, and 5A Ic = 0.003. 

^-correction to be subtracted from the ground-based magnitude to transform to the space-borne equivalent passband (see 
Tabled]). Mean uncertainties (mag) are 5SC B = 0.033, 5SC B = 0.041, 5SC Rc = 0.053, and SSC Ic = 0.014. 



